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We develop a method for the multifractal characterization of nonstationary time series, which is based on 
a generalization of the detrended fluctuation analysis (DFA). We relate our multifractal DFA method to the 
standard partition function-based multifractal formalism, and prove that both approaches are equivalent for 
stationary signals with compact support. By analyzing several examples we show that the new method can 
reliably determine the multifractal scaling behavior of time series. By comparing the multifractal DFA results 
for original series to those for shuffled series we can distinguish multifractality due to long-range correlations 
from multifractality due to a broad probability density function. We also compare our results with the wavelet 
transform modulus maxima (WTMM) method, and show that the results are equivalent. 
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I. INTRODUCTION 

In recent years the detrended fluctuation analysis (DFA) 
method [[!]] has become a widely-used technique for the de- 
termination of (mono-) fractal scaling properties and the de- 
tection of long-range correlations in noisy, nonstationary time 
series [||-^|]. It has successfully been applied to diverse fields 
such as DNA sequences heart rate dynamics [0^], neu- 
ron spiking [^[], human gait JlC|], long-time weather records 
[fill], cloud structure [O], geology [pj|], ethnology Jl4| ], eco- 
nomics time series []15|], and solid state physics JlM. One 
reason to employ the DFA method is to avoid spurious detec- 
tion of correlations that are artifacts of nonstationarities in the 
time series. 

Many records do not exhibit a simple monofractal scaling 
behavior, which can be accounted for by a single scaling ex- 
ponent. In some cases, there exist crossover (time-) scales 
Sx separating regimes with different scaling exponents [Dp]], 
e. g. long-range correlations on small scales s <C s x and an- 
other type of correlations or uncorrelated behavior on larger 
scales s 3> s x . In other cases, the scaling behavior is more 
complicated, and different scaling exponents are required for 
different parts of the series [||]. This occurs, e. g., when the 
scaling behavior in the first half of the series differs from the 
scaling behavior in the second half. In even more compli- 
cated cases, such different scaling behavior can be observed 
for many interwoven fractal subsets of the time series. In this 
case a multitude of scaling exponents is required for a full de- 
scription of the scaling behavior, and a multifractal analysis 
must be applied. 

In general, two different types of multifractality in time se- 
ries can be distinguished: (i) Multifractality due to a broad 
probability density function for the values of the time series. 
In this case the multifractality cannot be removed by shuf- 
fling the series, (ii) Multifractality due to different long-range 
(time-) correlations of the small and large fluctuations. In this 
case the probability density function of the values can be a 



regular distribution with finite moments, e. g. a Gaussian 
distribution. The corresponding shuffled series will exhibit 
non-multifractal scaling, since all long-range correlations are 
destroyed by the shuffling procedure. If both kinds of mul- 
tifractality are present, the shuffled series will show weaker 
multifractality than the original series. 

The simplest type of multifractal analysis is based upon the 
standard partition function multifractal formalism, which has 
been developed for the multifractal characterization of nor- 
malized, stationary measures fll7|-p0]]. Unfortunately, this 
standard formalism does not give correct results for nonsta- 
tionary time series that are affected by trends or that cannot 
be normalized. Thus, in the early 1990s an improved multi- 
fractal formalism has been developed, the wavelet transform 
modulus maxima (WTMM) method [|T]|, which is based on 
wavelet analysis and involves tracing the maxima lines in the 
continuous wavelet transform over all scales. Here, we pro- 
pose an alternative approach based on a generalization of the 
DFA method. This multifractal DFA (MF-DFA) does not re- 
quire the modulus maxima procedure, and hence does not in- 
volve more effort in programming than the conventional DFA. 

The paper is organized as follows: In Section II we de- 
scribe the MF-DFA method in detail and show that the scal- 
ing exponents determined via the MF-DFA method are iden- 
tical to those obtained by the standard multifractal formalism 
based on partition functions. In Section III we introduce sev- 
eral multifractal models, where the scaling exponents can be 
calculated exactly, and compare these analytical results with 
the numerical results obtained by MF-DFA. In Section IV, we 
show how the comparison of the MF-DFA results for original 
series to the MF-DFA results for shuffled series can be used to 
determine the type of multifractality in the series. In Section 
V, we compare the results of the MF-DFA with those obtained 
by the WTMM method for nonstationary series and discuss 
the performance of both methods for multifractal time series 
analysis. 
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II. MULTIFRACTAL DFA 

A. Description of the method 

The generalized multifractal DFA (MF-DFA) procedure 
consists of five steps. The first three steps are essentially iden- 
tical to the conventional DFA procedure (see e. g. Let 
us suppose that Xk is a series of length N, and that this series 
is of compact support, i.e. Xk — for an insignificant fraction 
of the values only. 

• Step 1: Determine the "profile" 

i 

Y(i) = J2lxk~(x)}, i = l,...,N. (1) 

k=i 

Subtraction of the mean (x) is not compulsory, since it would 
be eliminated by the later detrending in the third step. 

• Step 2: Divide the profile Y(i) into N s = int (iV/s) non- 
overlapping segments of equal length s. Since the length N of 
the series is often not a multiple of the considered time scale s, 
a short part at the end of the profile may remain. In order not to 
disregard this part of the series, the same procedure is repeated 
starting from the opposite end. Thereby, 2N S segments are 
obtained altogether. 

• Step 3: Calculate the local trend for each of the 2N S seg- 
ments by a least-square fit of the series. Then determine the 
variance 

F 2 (s, v) = -J2 Pit? - + i] - V^)} 2 (2) 
s L — ' 

»=i 

for each segment v, v = 1 , . . . , N a and 

1 s 

F 2 {s, v) = -Y\ {Y[N - {v - N s )s +t] - vA*)} 2 (3) 

for v = N s + 1, . . . , 2N S . Here, y v (i) is the fitting polynomial 
in segment v. Linear, quadratic, cubic, or higher order poly- 
nomials can be used in the fitting procedure (conventionally 
called DFA1, DFA2, DFA3, . . .) |JTJ|]. Since the detrending 
of the time series is done by the subtraction of the polynomial 
fits from the profile, different order DFA differ in their capa- 
bility of eliminating trends in the series. In (MF-)DFAm [mfh 
order (MF-)DFA] trends of order m in the profile (or, equiv- 
alently, of order m — 1 in the original series) are eliminated. 
Thus a comparison of the results for different orders of DFA 
allows one to estimate the type of the polynomial trend in the 
time series [f3jQ]. 

• Step 4: Average over all segments to obtain the qth order 
fluctuation function 

where, in general, the index variable q can take any real value 
except zero [p2|. For q = 2, the standard DFA procedure is 



retrieved. We are interested in how the generalized q depen- 
dent fluctuation functions F q (s) depend on the time scale s 
for different values of q. Hence, we must repeat steps 2 to 4 
for several time scales s. It is apparent that F q (s) will increase 
with increasing s. Of course, F q (s) depends on the DFA order 
m. By construction, F q (s) is only defined for s > rn + 2. 
• Step 5: Determine the scaling behavior of the fluctuation 
functions by analyzing log-log plots F q (s) versus s for each 
value of q. Several examples of this procedure will be shown 
in Section III. If the series Xi are long-range power-law corre- 
lated, F q (s) increases, for large values of s, as a power-law, 

F q (a)~a h M. (5) 

In general, the exponent h(q) may depend on q. For stationary 
time series, h(2) is identical to the well-known Hurst expo- 
nent H (see e. g. Jl7|]). Thus, we will call the function h(q) 
generalized Hurst exponent. 

For monofractal time series with compact support, h(q) is 
independent of q, since the scaling behavior of the variances 
F 2 (s, v) is identical for all segments v, and the averaging pro- 
cedure in Eq. (Q) will give just this identical scaling behavior 
for all values of q. Only if small and large fluctuations scale 
differently, there will be a significant dependence of h(q) on q: 
If we consider positive values of q, the segments v with large 
variance F 2 {v) (i. e. large deviations from the corresponding 
fit) will dominate the average F q {s). Thus, for positive values 
of q, h(q) describes the scaling behavior of the segments with 
large fluctuations. Usually the large fluctuations are character- 
ized by a smaller scaling exponent h(q) for multifractal series 
[p3|]. On the contrary, for negative values of q, the segments 
v with small variance F 2 (y) will dominate the average F q (s). 
Hence, for negative values of q, h(q) describes the scaling 
behavior of the segments with small fluctuations, which are 
usually characterized by a larger scaling exponent. 

However, the MF-DFA method can only determine posi- 
tive generalized Hurst exponents h(q), and it already becomes 
inaccurate for strongly anti-correlated signals when h(q) is 
close to zero. In such cases, a modified (MF-)DFA technique 
has to be used. The most simple way to analyze such data 
is to integrate the time series before the MF-DFA procedure. 
Hence, we replace the single summation in Eq. (|l|), which is 
describing the profile from the original data Xk, by a double 
summation, 

i 

f(i) = £[y(fc)-<F)]. (6) 

Following the MF-DFA procedure as described above, we ob- 
tain generalized fluctuation functions F q (s) described by a 
scaling law as in Eq. (|j), but with larger exponents h(q) — 
h(q) + l, 

F q (s) =s h ^ +1 . (7) 

Thus, the scaling behavior can be accurately determined even 
for h(q) which are smaller than zero (but larger than —1) for 
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some values of q. We note that F q (s)/s corresponds to F q (s) 
in Eq. If we do not subtract the average values in each step 
of the summation in Eq. (^|), this summation leads to quadratic 
trends in the profile Y(i). In this case we must employ at least 
the second order MF-DFA to eliminate these artificial trends. 



B. Relation to standard multifractal analysis 

For stationary, normalized records with compact support 
the multifractal scaling exponents h(q) defined in Eq. (||) are 
directly related, as shown below, to the scaling exponents r(q) 
defined by the standard partition function-based multifractal 
formalism. 

Suppose that the series Xk of length N is a stationary, nor- 
malized sequence. Then the detrending procedure in step 3 
of the MF-DFA method is not required, since no trend has to 
be eliminated. Thus, the DFA can be replaced by the stan- 
dard Fluctuation Analysis (FA), which is identical to the DFA 
except for a simplified definition of the variance for each seg- 
ment v, v = 1, . . . , N s , in step 3 [see Eq. (0)]: 



F* A (s,v) = [Y(vs)-Y((v-l)s)f 



(8) 



Inserting this simplified definition into Eq. (^J) and using 
Eq. (0), we obtain 



(9) 



For simplicity we can assume that the length N of the series 
is an integer multiple of the scale s, obtaining N s = N/ s and 
therefore 



N/s 



^W-y(( v -i»|« 



(10) 



This already corresponds to the multifractal formalism used 
e. g. in [ p"8"p0| ] . In fact, a hierarchy of exponents H q similar 
to our h{q) has been introduced based on Eq. ( |l(j[ ) in [[T|]. 

In order to relate also to the standard textbook box counting 
formalism flT^ , p"9[ l, we employ the definition of the profile in 
Eq. dl|). It is evident that the term Y(vs) - Y((v - l)s) in 
Eq. (HQ) is identical to the sum of the numbers xu within each 
segment v of size s. This sum is known as the box probability 
p s {v) in the standard multifractal formalism for normalized 
series Xk, 



VS 

£ 

fc=(i/-l)s+l 



x k = Y{vs)-Y{{v-l)s). (11) 



The scaling exponent r(q) is usually defined via the partition 
function Z q (s), 



N/s 



Z q (s)=J2\Ps(v)\ q ~J 



(?) 



(12) 



where q is a real parameter as in the MF-DFA above. Some- 
times r(q) is defined with opposite sign (see e. g. JT^]). 

Using Eq. ( pi] ) we see that Eq. ( p"2| ) is identical to Eq. (|To|), 
and obtain analytically the relation between the two sets of 
multifractal scaling exponents, 



r(q) = qh(q) - 1. 



(13) 



Thus, we have shown that h(q) defined in Eq. (||) for the MF- 
DFA is directly related to the classical multifractal scaling ex- 
ponents r{q). Note that h(q) is different from the generalized 
multifractal dimensions 



D(q) 



r(q) qh(q) - 1 



9-1 



g-i 



(14) 



that are used instead of r(q) in some papers. While h(q) is 
independent of q for a monofractal time series with compact 
support, D(q) depends on q in that case. Our assumption of 
compact support of the series Xk can be directly observed in 
Eq. (pi]), since the fractal dimension of the support is D(0) = 
-r(0) = l. 

Another way to characterize a multifractal series is the sin- 
gularity spectrum f(a), that is related to r(q) via a Legendre 
transform 



r'(q) and f(a)=qa-T(q). 



(15) 



Here, a is the singularity strength or Holder exponent, while 
/(a) denotes the dimension of the subset of the series that is 
characterized by a. Using Eq. (pj|), we can directly relate a 

and f(a) to h(q), 

a = h(q) + qh'{q) and f{a) = q[a - h(q)] + 1. (16) 



in. FOUR ILLUSTRATIVE EXAMPLES 

A. Example 1: monofractal uncorrelated and long-range 
correlated series 

As a first example we apply the MF-DFA method to 
monofractal series with compact support, for which the gen- 
eralized Hurst exponent h(q) is expected to be independent of 

q, 



h(q) = H and r(q) = qH - 1. 



(17) 



Such series have been discussed in the context of conventional 
DFA in several studies before, see e.g. [pj-pl]- Long-range cor- 
related random numbers are usually generated by the Fourier 
transform method, see e. g. jl7| , p4[ |. Using this method we can 
generate long-range anti-correlated (0 < H < 0.5), uncorre- 
lated (H = 0.5), or (positively) long-range correlated (0.5 < 
H < 1) series. The latter are characterized by a power-law de- 
cay of the autocorrelation function C(s) = (x^ Xk+ S ) ~ s~ 7 
for large scales s with 7 = 2 — 2H if the series is stationary. 
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Alternatively, all stationary long-range correlated series can 
be characterized by the power-law decay of their power spec- 
tra, S(f) ~ f~ f3 with frequency / and /3 = 2H — 1. Note 
that H corresponds to the Hurst exponent of the integrated 
time series here. 

Figure [j] shows the generalized fluctuation functions F q (s) 
for all three types of monofractal series (H = 0.75, 0.5, 0.25) 
and several q values. On large scales s, we observe the 
expected power-law scaling behavior according to Eq. (^), 
which corresponds to straight lines in the log-log plot. In 
Fig. |l](d), the scaling exponents h(q) determined from the 
slopes of these straight lines are shown versus q. Although 
a slight q dependence is observable, the values of h(q) are al- 
ways very close to the H of the generated series that has been 
analyzed. The degree of the q dependence observed for this 
monofractal series allows to estimate the usual fluctuation of 
h(q) to be expected for monofractal series in general. 

Next, we analyze multifractal series for which r(q) can be 
calculated exactly, and compare the numerical results with the 
expected scaling behavior. 



B. Example 2: binomial multifractal series 

In the binomial multifractal model [jj^UjlJ, a series of N 
2"max num bers k with k = 1, . . . , N is defined by 



= "(*-l)( 1 _ )nm«-n(*-l) j 



(18) 



where 0.5 < a < 1 is a parameter and n(k) is the number of 
digits equal to 1 in the binary representation of the index k, 
e. g. n(13) = 3, since 13 corresponds to binary 1101. 

The scaling exponents r(q) can be calculated straightfor- 
wardly. According to Eqs. ( p"l| ) and (|l8]) the box probability 
P2s{v) in the vth segment of size 2s is given by 

P2s{v) = Ps(2v - 1) + Ps{2v) 

= [(1 - a)/a + l]p s {2v) = p s {2v)/a. 
Thus, according to Eqs. ( |l2| ) and (|l8|), 

N/s N/2s 
Z q { S ) = £>.(!/)]« = £ iPsi^ I)]' + iPs(^)} q 



ai 



-. N/2s 



1 



E ^(^) 



N/2s 



[(1 - aY + a"] J2 [Pis{v)\ q = [(1 - «) 9 + a 9 ] Z q {2s) 



and according to Eqs. ( |l2| ) and (pj|), 

■ln{a q + (1 - o)«] 



r(«) 



ln(2) 

1 _ ln[a 9 + (1 - o)«] 
9 



gln(2) 



(19) 
(20) 



Note that r(0) = — 1 as required. There is a strong non- 
linear dependence of r(q) upon q, indicating multifractality. 
The same information is comprised in the q dependence of 
h(q). The asymptotic values are h(q) — > — ln(a)/ln(2) for 
q — » +c» and h(q) — > — ln(l — a)/ ln(2) for q — > — oo. They 
correspond to the scaling behavior of the largest and weakest 
fluctuations, respectively. Note that h(q) becomes indepen- 
dent of q in the asymptotic limit, while r(q) approaches linear 
q dependences. 

Figure ^ shows the MF-DFA fluctuation functions F q (s) 
for the binomial multifractal model with a — 0.75. The re- 
sults for MF-DFA 1 and MF-DFA4 are compared in parts (a) 
and (b). Fig. ||(c) shows the corresponding slopes h(q) for 
three values of a together with the exact results obtained from 
Eq. (M). The numerical results are in good agreement with 
Eq. (|20[), showing that the MF-DFA correctly detects the mul- 
tifractal scaling exponents. Figures ^(d) and (e) show the cor- 
responding exponents r(q) = qh(q) — 1 [see Eq. (|l3"|)] and 
the corresponding f(a) spectrum calculated from h(q) using 
the modified Legendre transform (|l6|). Both are also in good 
agreement with Eq. (|l9|). We have also checked that the re- 
sults for the binomial multifractal model remain unchanged if 
the double summation technique [see Eq. (^)] is applied. We 
obtain slopes h(q) = h(q) + 1 as expected in Eq. (Q). Note 
that there is no need to use this modification, except if h(q) is 
close to zero or has negative values. 



C. Example 3: dyadic random cascade model with log-Poisson 
distribution 



For another independent test of the MF-DFA, we employ an 
algorithm based on random cascades on wavelet dyadic trees 
proposed in p5| ] (see also p6|]). This algorithm builds a ran- 
dom multifractal series by specifying its discrete wavelet co- 



efficients c nm , defined recursively, 



ci.i = 1, 



C'Ti , 2 rn 



where n — 2,...,n max (with = 2™ max ) and m 

1 2™~ 2 . The values of W are taken from a log-Poisson 

distribution, \ W\ = exp(P In S + 7), where P is Poisson dis- 
tributed with (P) = A. There are three independent parame- 
ters, A, 5, and 7. Inverse wavelet transform is applied to create 
the multifractal random series Xk once the wavelet coefficients 
c n m are known, 



EE' 

n—l m=l 



(21) 



where ip n . m (k) is a set of wavelets forming an orthonor- 
mal wavelet basis. Here, we employ the Haar wavelets, 
ip n , m (k) = 2( n " n >»»*- 1 )/ 2 V>[2 n - n "»~ - x fc - to] with ip(x) = 
lforO < x < 0.5,ip(x) = -1 for 0.5 < x < land tp(x) = 
otherwise. For this model the multifractal scaling exponents 
are given by [S 
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h(q) = [\(l-5«)- 7 q}/(qhi2). 



(22) 
(23) 



Figure || shows the MF-DFA fluctuation functions F q (s) 
for the dyadic random cascade model. The numerically deter- 
mined slopes h(q) for three sets of parameters are compared 
with the exact results obtained from Eq. (Q) and the good 
agreement shows that the MF-DFA correctly detects the mul- 
tifractal scaling exponents. Large deviations occur only for 
very small moments (q < —10), indicating that the range of q 
values should not exceed —10. 



D. Example 4: uncorrelated multifractal series with power-law 
distribution function 



The examples discussed in the previous three subsections 
were based on series involving long-range correlations. In the 
present example we want to apply the MF-DFA method to 
an uncorrelated series, that nevertheless exhibits multifractal 
scaling behavior due to the broad distribution of its values. We 
denote by P(x) the probability density function of the values 
Xk in the series. The distribution P(x) does not affect the 
multifractality of a series on large scales s, if all moments 



\x\ q P{x)dx 



(24) 



are finite. Here we choose a (normalized) power-law proba- 
bility distribution function, 



P(x) 



ax 



" (Q+1) for 1 < x < oo with a > (25) 



and P(x) = for x < 1, where already the second moment 
diverges if a < 2. In this case, the series exhibits multifractal 
scaling behavior on all scales. Note, that Eq. (|5j) becomes 
identical to a Levy distribution of class a for large values of 
x. The parameter a is not related to the Holder exponent a in 
Eq. (p3|). The scaling properties of random walks with Levy 
distributed steps (Levy flights and Levy walks) have been an- 
alyzed in [p7[-p9|]. The multifractal nature of Levy processes 
has been investigated in [ pOpip . 

In order to derive the multifractal spectrum, let us consider 
s uncorrelated random numbers 77., k = 1, . . . , s, distributed 
homogeneously in the interval [0, 1]. Obviously, the typical 
value of the minimum of the numbers, r m j n (s) = mm s k=1 r k , 
will be r m i n (s) = 1/s. It can be easily shown that the num- 
bers r k are transformed into numbers Xk distributed accord- 
ing to the power-law probability distribution function ( ^sj ) by 
Tk — * Xk = r~ k 1 . Thus, the typical value of the maximum 
of the x k will be x max (s) = ma,x s k=1 x k = [r lriin {s)}- 1/a = 

If a < 2, the fluctuations of the profile Y[i) [Eq. (jl|)] and 
the corresponding DFA variance F 2 {s,v) [Eq. (|)] will be 
dominated by the square of the largest value i^fs) = s 2 l a 
in the segment of s numbers, since the second moment of the 



distribution ( |25[ ) diverges. Now the whole series consists of 
N s = int(7V/s) segments of length s and not just of one 
segment. For some segments v, [F 2 {s,v)] 1 / 2 is larger than 
its typical value x ma ^(s) — s 1 /™, since the maximum within 
the whole series of length N is x ma ^(N) = N 1 / . In or- 
der to calculate F q (s) [Eq. (Q)], we need to take into account 
the whole distribution P s (y) of the values y = [F 2 (s, v)} 1 ^ 2 . 
Since each of the maxima in the N s segments corresponds 
to an actual number x k and these x k are random numbers 
from the power-law distribution (^), it becomes obvious, that 
the distribution of the maxima will have the same form, i. e. 
Ps{y) ~ P(x = y) for large y. Small values of y are ex- 
cluded because of the maximum procedure, but the large x k 
values are very likely to be identical to the maxima of the 
corresponding segments. Since the smallest maxima for seg- 
ments of length s are of the order of x max (s) = s 1 /", the 
lower cutoff for P s (y) must be proportional to s 1 /". From 
the normalization condition J A 1/a P s (y) dy = 1 (with an 
unimportant prefactor A < 1) we get 



P s (y) = A a asy^ a+1 l 



(26) 



Now F q (s) [Eq. (Q)] can be calculated by integration from 
the minimum value As x l a of y = [F 2 (s, v)] 1 / 2 to the maxi- 
mum value N 1 / . For s N we obtain 



F g (s) 



y q P s {y) dy 



1/9 



- A a sN q l a - 1 - A q s q/a 

Comparing with Eq. (Q), we finally get 

h(a)~l 1/q ( q>a) 
n[Q) Ua (q<a) 



As 1 /" 

1 ' q ( s^ q (q>a) 
s 1 /" (q < a) 



(27) 



Note that r(q) follows a linear q dependence, r(q) — q/a — 1 
for q < a, while it is equal to zero for q > a according 
to Eq. dl3[). Hence, the series of uncorrelated power-law dis- 
tributed values has rather bi-fractal [|TJ| instead of multifractal 
properties. Since h(2) = 1/2 holds exactly for all values of 
a, it is not possible to recognize the multifractality due to the 
broad power-law distribution of the values if only the conven- 
tional DFA is applied. The second moment shows just the 
uncorrelated behavior of the values. In a very recent preprint 
J^ ] this behavior has been interpreted as a failure of the DFA 
and corresponding non-detrending methods for series with a 
broad distribution, and another method to determine the ex- 
ponent 1 1 a has been proposed. We believe that a multifractal 
description with more than one exponent is required to charac- 
terize this kind of series, and thus any method calculating just 
one exponent will be insufficient for a full characterization. 

Figure |](a) shows the MF-DFA3 fluctuation functions for 
series of independent random numbers x k G [1, 00) dis- 
tributed according to Eq. (25) with a = 1. Since the scaling 
exponents h(q) become very close to zero asymptotically for 
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large positive values of q according to Eq. (|27|), we must use 
the modified MF-DFA technique involving the double sum as 
described in the last paragraph of Subsection II. A. Hence, for 
this technical reason, F q (s)/s is calculated instead of F q (s). 
The corresponding slopes h(q) — 1 are identical to h(q), see 
Eq. (0). In Fig. 0(b) the slopes h{q) for series with a = 0.5, 
1.0, and 2.0 are compared with the theoretical result Eq. (j27|), 
and nice agreement is observed. 

IV. COMPARISON OF THE MULTIFRACTALITY FOR 
ORIGINAL AND SHUFFLED SERIES 

A. Distinguishing the two types of multifractality 

As already mentioned in the introduction, two different 
types of multifractality in time series can be distinguished. 
Both of them require a multitude of scaling exponents for 
small and large fluctuations, (i) Multifractality of a time se- 
ries can be due to a broad probability density function for the 
values of the time series, and (ii) multifractality can also be 
due to different long-range correlations for small and large 
fluctuations. The example discussed in Subsection III.D, the 
uncorrelated multifractal series with a power-law probability 
density function, is of type (i), while the examples discussed 
in Subsections III. A - III.C are of type (ii), where the prob- 
ability density function of the values is a regular distribution 
with finite moments [[32]] . 

Now, we would like to distinguish between these two types 
of multifractality. The most easy way to do so is by analyzing 
also the corresponding randomly shuffled series. In the shuf- 
fling procedure the values are put into random order, and thus 
all correlations are destroyed. Hence the shuffled series from 
multifractals of type (ii) will exhibit simple random behavior, 
ft'shuf (?) = 0.5, i. e. non-multifractal scaling like in Fig. 0(b). 
For multifractals of type (i), on the contrary, the original h(q) 
dependence is not changed, h(q) = h s \ lu f(q), since the mul- 
tifractality is due to the probability density, which is not af- 
fected by the shuffling procedure. If both kinds of multifrac- 
tality are present in a given series, the shuffled series will show 
weaker multifractality than the original one. 

The effect of the shuffling procedure is illustrated in 
Fig. 0(a), where the MF-DFA2 fluctuation functions F!% { (s) 
and Ff^ (s) are shown for shuffled series for three of the 
multifractal examples taken from the previous section. Ran- 
dom behavior, /i s huf (<?) = 0.5, is observed for the series that 
were long-range correlated or generated from the dyadic ran- 
dom cascade model before the shuffling procedure [upper four 
curves in Fig. 0(a)]. In contrast, we observe the original multi- 
fractal scaling for the shuffled multifractal series with power- 
law probability density function P{x) ~ x~ 2 [lower two 
curves in Fig. 0(a)]. The /i s huf (<?) dependences are shown in 
Fig. 0, which can be compared with the corresponding slopes 
shown in Figs. 0(d), |(b), and 0(b). Thus, the fluctuation 
analysis of the shuffled series, F* (s), directly indicates the 



presence of type (i) multifractality, which is due to a broad 
probability distribution, by deviations from /i 8 h u f (?) = 0.5. 

Now we want to determine directly the magnitude of the (ii) 
multifractality, which is due to correlations. For that purpose 
we compare the fluctuation function for the original series, 
F q (s), with the result for the corresponding shuffled series, 
F q huf (s). Differences between these two fluctuation func- 
tions directly indicate the presence of correlations in the orig- 
inal series. These differences can be observed best in a plot 
of the ratio F q (s) / F q hu! (s) versus s jf33|]. Since the anoma- 
lous scaling due to a broad probability density affects F q (s) 
andF| huf (s) in the same way, only multifractality due to cor- 
relations will be observed in F q (s)/F^ (s). This is illus- 
trated in Fig. 0(b) for the same three multifractal examples as 
in Fig. 0Ja). In order not to have increased statistical errors in 
the results when considering the ratio F q (s) / F q huf (s) instead 
of F q (s) itself, F q (s) can be calculated by averaging over 
a large number of randomly shuffled series generated from the 
same original series. 

The scaling behavior of the ratio is 

Fq(s)/F q huf ( s ) ~ s h ^~ h « h ^ = S ' W ^. (28) 

Note that h(q) — /i s huf(<?) + h COI (q). If only distribution 
multifractality [type (i)] is present, h(q) — /i s huf (<?) depends 
on q and h COI (q) — 0. On the other hand, deviations of 
hcoi(q) from zero indicate the presence of correlations, and 
a q dependence of h cm (q) indicates correlation multifractal- 
ity [type (ii)]. If only correlation multifractality is present, 
h s hui(q) — 0.5 and h(q) = 0.5 + h COI (q). If both, distribu- 
tion multifractality and correlation multifractality are present, 
both, ftshuf (q) an d h cm (q) depend on q. 

B. Significance of the results 

In Figs. 0^ we have shown the results of the MF-DFA for 
single configurations of long time series. Now we address 
the significance and accuracy of the MF-DFA results for short 
series. How much do the numerically determined exponents 
h(q) vary from one configuration (sample series) to the next, 
and how close are the average values to the theoretical val- 
ues? In other words, how large are the statistical and system- 
atical deviations of exponents practically determined by the 
MF-DFA for finite series? These questions are particularly 
important for short series, where the statistics is poor. If the 
values of h(q) are determined inaccurately, the multifractal 
properties will be reported inaccurately or even false conclu- 
sions on multifractal behavior might be drawn for monofractal 
series. 

To address the significance and accuracy of the MF-DFA 
results we generate, for each of the three examples considered 
already in Fig. 100 series of length TV = 2 13 = 8192 and 
calculate h(-10), h(+10), h shul (~10), and /i slm f(+10) for 
each of these series. The corresponding histograms are shown 
in Fig. [7| For the long-range power-law correlated series with 
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H = 0.75 we find the following mean values and standard 
deviations of the generalized Hurst exponents: 

/i(-10) = 0.80 ± 0.03, /i s huf(-10) = 0.56 ± 0.02, 
/i(+10) = 0.72 ± 0.04, h shui (+10) = 0.48 ± 0.02. 

The mean values for the original series are rather close to, but 
not identical to the theoretical value H = 0.75. The mean 
value for q = —10 is about two standard deviations larger 
than 0.75, while the value for q = +10 is slightly smaller. 
These deviations, though, certainly cannot indicate multifrac- 
tality, since we analyzed monofractal series. Instead, they are 
due to the finite, random series, where parts of the series have 
slightly larger and slightly smaller scaling exponent just by 
statistical fluctuations. By considering negative values of q 
we focus on the parts with small fluctuations, which are usu- 
ally described by a larger scaling exponent [^3[]. For positive 
values of q we focus on the parts with large fluctuations usu- 
ally described by a smaller value of h. Thus for short records 
we always expect a slight difference between h(— 10) and 
h(+10) even if the series are monofractal. If this difference 
is weak, one has to be very careful with conclusions about 
multifractality. Practically it is always wise to compare with 
generated monofractal series with otherwise similar properties 
before drawing conclusions regarding the multifractality of a 
time series. In addition to the statistical fluctuations of the h 
values, the average h(— 10) is usually determined slightly too 
large, while h(+10) is slightly too small. The same behav- 
ior is obtained if the WTMM method is used instead of the 
MF-DFA, as we will show in Subsection V.C. 

The same kind of difference is also observed for the aver- 
age /ighuf (— 10) and ft. s huf (+10) values. After all correlations 
have been destroyed by the shuffling, /i snu f = 0.5 is expected 
since the probability density is Gaussian with all finite mo- 
ments. The deviations from /i s huf = 0.5 we observe for the 
finite random series are characteristic for monofractal series 
of this length (N = 8192). Only for the second moment we 
obtain /i s h u f (2) = 0.5 exactly if a sufficient number of series 
is considered. 

For multifractal series generated from the dyadic random 
cascade model, Fig. 0(c,d) shows the histograms of the scal- 
ing exponents h(— 10), h(+l0), /i s huf (— 10), and /ishuf (+10). 
Their averages and standard deviations, 

/i(-10) = 0.69 ± 0.04, /i shu f (-10) = 0.57 ± 0.02, 
h(+10) = 0.54 ± 0.02, /i s huf (+10) = 0.48 ± 0.02. 

have to be compared with the theoretical values from Eq. (|23|), 
/i(-10) = 0.743 and h(+W) = 0.567. Surprisingly, the 
mean h(— 10) is smaller than the theoretical value in this ex- 
ample, but for the mean h(+10) the deviation is similar to 
the deviation observed for the monofractal data in the pre- 
vious example. Again, similar results are obtained with the 
WTMM method. For the shuffled series, the mean general- 
ized Hurst exponents are practically identical to those for the 
shuffled monofractal series [the average /i s huf (— 10) is larger 



by half the standard deviation], and both are evidently consis- 
tent with monofractal uncorrected behavior, h(q) — 0.5, as 
discussed above. Hence, the series from the dyadic random 
cascade model show no signs of distribution multifractality 
and are characterized by correlation multifractality only. 

The histograms of the scaling exponents for our last exam- 
ple, the power-law distributed random numbers with P(x) ~ 
x~ 2 , are shown in Fig. [7](e,f). The corresponding mean values 
and standard deviations, 

/i(-10) = 1.24 ± 0.09, /i shu f (-10) = 1.26 ± 0.09, 
/i(+10) = 0.11 ± 0.03, ft^ (+10) = 0.11 ± 0.04. 

show obviously no differences between original and shuffled 
series as expected for uncorrelated series. This indicates that 
the multifractality is due to the broad probability density func- 
tion only. The values have to be compared with h(— 10) = 1 
and /i(+10) = 0.1 from Eq. (^7|). As usual, the average value 
of h(— 10) is too large because we analyzed short series. 

V. COMPARISON TO THE WAVELET TRANSFORM 
MODULUS MAXIMA METHOD 

A. Brief description of the wavelet transform modulus maxima 
method 

The wavelet transform modulus maxima (WTMM) method 
J2l| ] is a well-known method to investigate the multifractal 
scaling properties of fractal and self-affine objects in the pres- 
ence of nonstationarities. It is an application of the wavelet 
transform with continuous basis functions. One defines the 
wavelet-transform of a series x k of length N by 

1 N 

W{n,s) = -Y^x k ip[{k-n)/s], (29) 
s ^— ' 
fe=i 

Note that in this case the series x k are analyzed directly in- 
stead of the profile Y(i) defined in Eq. (jl|). Here, the function 
ip{x) is the analyzing wavelet and s is, as above, the scale 
parameter. The wavelet is chosen orthogonal to the possi- 
ble trend. If the trend can be represented by a polynomial, 
a good choice for ijj(x) is the m-th derivative of a Gaussian, 
-0( m )(x) = d m (e~ x2/2 )/dx m . This way, the transform elimi- 
nates trends up to (in — l)th order. 

Now, instead of averaging over all values of W(n, s), one 
averages, within the modulo-maxima method, only the lo- 
cal maxima of |T^(n, s)|. First, one determines for a given 
scale s, the positions rt,; of the local maxima of |VK(n, s)| 
as function of n, so that |W(nj — l,s)\ < \W(rii,s)\ > 
\W(rii + 1, s)\ for i = 1, . . . , i m ax- Then one sums up the 
qth power of these maxima, 

i max 

Z{q,s) = g |W(ni,s)|«. (30) 
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The reason for this maxima procedure is that the absolute 
wavelet coefficients |W(n, s)| can become arbitrarily small. 
The analyzing wavelet ip(x) must always have positive values 
for some x and negative values for other x, since it has to be 
orthogonal to possible constant trends. Hence there are always 
positive and negative terms in the sum (p9|), and these terms 
might cancel. If that happens, |VK(n, s)\ can become close 
to zero. Since such small terms would spoil the calculation 
of negative moments in Eq. ([30|), they have to be eliminated 
by the maxima procedure. In the MF-DFA, the calculation of 
the variances F 2 (s, v) in Eq. (||), i. e. the deviations from the 
fits, involves only positive terms under the summation. The 
variances cannot become arbitrarily small, and hence no max- 
imum procedure is required for series with compact support. 

In addition, the MF-DFA variances will always increase if 
the segment length s is increased, because the fit will always 
be worse for a longer segment. In the WTMM method, in 
contrast, the absolute wavelet coefficients s) \ need not 

increase with increasing scale s, even if only the local maxima 
are considered. The values ^(rt, s)\ might become smaller 
for increasing s since just more (positive and negative) terms 
are included in the summation (|9|), and these might cancel 
even better. Thus, an additional supremum procedure has been 
introduced in the WTMM method in order to keep the depen- 
dence of Z(q, s) on s monotonous: If, for a given scale s, a 
maximum at a certain position rii happens to be smaller than a 
maximum at r n! i w rii for a lower scale s' < s, then W(rii, s) 



is replaced by W^, s') in Eq. (|30|). There is no need for such 
a supremum procedure in the MF-DFA. 

Often, scaling behavior is observed for Z(q, s), and scal- 
ing exponents r(q) can be defined that describe how Z(q, s) 
scales with s, 



Z(q,s) 



b t(q) 



(31) 



The exponents r(q) characterize the multifractal properties of 
the series under investigation, and theoretically they are iden- 
tical to the r(q) defined in Eq. ( |l2| ) and related to h(q) in 
Eq. ([l3]). 



B. Examples for series with nonstationarities 

Since the WTMM method has been developed to analyze 
multifractal series with nonstationarities, such as trends or 
spikes, we will compare its performance with the performance 
of the MF-DFA for such nonstationary series. In Fig. |] the 
MF-DFA fluctuation function F q (s) and its scaling behavior 
are compared with the rescaled WTMM partition sum Z(q, s) 
for the binomial multifractal described in Subsection III.B. 
To test the detrending capability of both methods, we have 
added linear as well as quadratic trends to the generated mul- 
tifractal series. The trends are removed by both methods, if 
a sufficiently high order of detrending is employed. The de- 
viations from the theoretical values of the scaling exponents 
h(q) [given by Eq. (Eoh] are of similar size for the MF-DFA 



and the WTMM method. Thus, the detrending capability and 
the accuracy of both methods is equivalent. 

We also obtain similar results for a monofractal long-range 
correlated series with additional spikes (outliers) that consist 
of large random numbers and replace a small fraction of the 
original series in randomly chosen positions. The spikes lead 
to multifractality on small scales s, while the series remains 
monofractal on large scales. Thus, the effects of the spikes 
are eliminated neither by the WTMM method nor by the MF- 
DFA, but both methods again give rather equivalent results. 



C. Significance of the results 

The last problem we address is a comparison of the signifi- 
cance of the results obtained by the MF-DFA and the WTMM 
method. The significance of the MF-DFA results has already 
been discussed in detail in Subsection IV.B. Here we will 
compare the significance of both methods for short and long 
series. 

We begin with the significance of the results for random 
series involving neither correlations nor a broad distribution 
[as in Fig. [lib)]. Fig. ^ shows the distribution of the mul- 
tifractal Hurst exponents h(— 10) and h(+10) calculated by 
the MF-DFA as well as by the WTMM using the relation 
M?) = [t(<1) + !]/<? based on Eq. (|l3|). Similar to the results 
presented in Fig. |j, we have analyzed 100 generated series of 
uncorrelated random numbers. In addition, we compare the 
results for the (relatively short) series length N — 2 13 = 8192 
and for N = 2 16 = 65532. Ideally, both, /i(-10) and 
h(+10), should be equal to the Hurst exponent of the uncorre- 
lated monofractal series, H = 0.5. The histograms show that 
similar deviations as well as remarkable fluctuations of the 
exponents occur for both methods, as discussed in Subsection 
IV.B for the MF-DFA. We find the following mean values and 
standard deviations, 



ft(-10) 



and h(+10) 



0.55 ±0.03 
0.52 ±0.02 
0.58 ±0.05 
0.56 ±0.03 
0.49 ±0.03 
0.49 ±0.02 
0.46 ±0.04 
0.48 ±0.02 



for MF-DFA (N = 8k) 
for MF-DFA (iV = 64/fc) 
for WTMM (N = 8k) 
for WTMM (N = 64fc) 
for MF-DFA (N = 8k) 
for MF-DFA (N = 64/fc) 
for WTMM (N = 8k) 
for WTMM (JV = 64fc) 



As already discussed in Subsection IV.B, the deviations of the 
average h(q) values from H — 0.5 do not indicate multifrac- 
tality. For the WTMM method and short series, one has to be 
very careful in order not to draw false conclusions from results 
like /i(-10) = 0.58 and fo(+10) = 0.46. The corresponding 
results of the MF-DFA are closer to the theoretical value. 

Figure [To] shows the distribution of the multifractal scaling 
exponents h(—10) and h(+10) calculated for generated mul- 
tifractal series from the binomial model with a = 0.75 de- 
scribed in Subsection III.B. Like for Fig. ^ 100 generated se- 
ries have been analyzed for each of the histograms. Now, the 
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differences between the distributions of h(— 10) and h(+10) 
are much larger, indicating multifractality. We find 



h(-W) 



and h(+10) 



1.88 ±0.06 

1.89 ±0.03 
1.86 ±0.05 
1.89 ±0.02 

0.50 ±0.02 
0.51 ±0.01 
0.46 ±0.01 
0.47 ±0.01 



for MF-DFA (N = 8k) 
for MF-DFA (N = 64k) 
for WTMM (N = 8k) 
for WTMM (TV = 64fc) 
for MF-DFA (N = 8k) 
for MF-DFA (N = 64k) 
for WTMM (N = 8k) 
for WTMM (TV = 64k) 



These values must be compared with the theoretical values 
h(-10) = 1.90 and ft,(+10) = 0.515 from Eq. @. Again, 
the MF-DFA results turn out to be slightly more significant 
than the WTMM results. The MF-DFA seems to have slight 
advantages for negative q values and short series, but in the 
other cases the results of the two methods are rather equiva- 
lent. Besides that, the main advantage of the MF-DFA method 
compared with the WTMM method lies in the simplicity of 
the MF-DFA method. 



Subsection V.A). There is no need, though, to employ a con- 
tinuously sliding window or to calculate the supremum over 
all lower scales for the MF-DFA, since the variances F 2 (s,v), 
which are determined by the deviations from the fit, will al- 
ways increase when the segment size s is increased. In the 
maxima MF-DFA procedure the generalized Hurst exponent 
h(q) defined in Eq. (Q) will depend on q and even diverge 
for q — > for monofractal series with non-compact support. 
Thus, it is more appropriate to consider the scaling exponent 
r(q), calculating 

[F 2 {s.u)] q l 2 ~ s T( - q \ (32) 

F 2 (s,v-l)<F 2 (s,v)>F 2 (s,y+l) 

This extended MF-DFA procedure will also be applicable for 
data with fractal support. 

In a later work we will apply the MF-DFA method to a 
range of physiological and meteorological data. 
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VI. CONCLUSION 

We have generalized the DFA, widely recognized as a 
method to analyze the (mono-) fractal scaling properties of 
nonstationary time series. The MF-DFA method allows a re- 
liable multifractal characterization of multifractal nonstation- 
ary time series. The implementation of the new method is not 
more difficult than that of the conventional DFA, since just 
one additional step, a q dependent averaging procedure, is re- 
quired. We have shown for stationary signals that the gen- 
eralized (multifractal) scaling exponent h(q) for series with 
compact support is directly related to the exponent r(q) of 
the standard partition function-based multifractal formalism. 
Further, we have shown in several examples that the MF-DFA 
method can reliably determine the multifractal scaling behav- 
ior of the time series, similar to the WTMM method which is 
a more complicated procedure for this purpose. For short se- 
ries and negative moments, the significance of the results for 
the MF-DFA seems to be slightly better than for the WTMM 
method. 

Contrary to the WTMM method, the MF-DFA method as 
described in Subsection II. A requires series of compact sup- 
port, because the averaging procedure in Eq. (Q) will only 
work if F 2 (s : v) > for all segments v. Although most 
time series will fulfill this prerequisite, it can be overcome by 
a modification of the MF-DFA technique in order to analyze 
data with fractal support: We restrict the sum in Eq. (Q) to the 
local maxima, i. e. to those terms F 2 (s, v) that are larger than 
the terms F 2 (s, v — 1) and F 2 (s, v + 1) for the neighboring 
segments. By this restriction all terms F 2 (s, v) that are zero 
or very close to zero will be disregarded, and series with frac- 
tal support can be analyzed. The procedure reminds slightly 
of the modulus maxima procedure in the WTMM method (see 
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FIG. 1. The MF-DFA fluctuation functions F q (s) are shown 
versus the scale s in log-log plots for (a) long-range correlated 
monofractal series with H — 0.75, (b) uncorrelated random series 
with H = 0.5 (white noise), and (c) long-range anti-correlated se- 
ries with H — 0.25. The different symbols correspond to the differ- 
ent values of the exponent q in the generalized averaging procedure, 
q = -10 (□), -2 (O), -0.2 (A), +0.2 (v), +2 (O), and q = +10 
(+). MF-DFA2 has been employed, and the curves have been shifted 
by multiple factors of 4 for clarity. The straight dashed lines have 
the corresponding slopes H and are shown for comparison. Part (d) 
shows the q dependence of the asymptotic scaling exponent h(q) de- 
termined by fits in the regime 200 < s < 5000 for H = 0.25 (□), 
0.5 (O), and 0.75 (A). The very weak dependence on q is consis- 
tent with monofractal scaling. In part (e) the q dependence of r(q), 
r(g) = qh(q) — 1, is shown. 



FIG. 2. The MF-DFA fluctuation functions F q (s) are shown ver- 
sus the scale s in log-log plots for the binomial multifractal model 
with a = 0.75 (a) for MF-DFA 1 and (b) for MF-DFA4. The symbols 
are the same as for Fig. 1 . The straight dashed lines have the corre- 
sponding theoretical slopes h(— 10) = 1.90 and h(+10) = 0.515 
and are shown for comparison. In part (c) the q dependence of the 
generalized Hurst exponent h(q) determined by fits in the regime 
50 < s < 500 is shown for MF-DFA 1 and a = 0.9 (A), a = 0.75 
(O), and a = 0.6 (□), as well as for MF-DFA4 and a = 0.75 (v). 
Parts (d) and (e) show the corresponding exponents r(q) and the cor- 
responding singularity spectrum f(a) for a = 0.75 determined by 
the modified Legendre transform (|lq), respectively. The lines are the 
theoretical values obtained from Eq. fl20J). 
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FIG. 3. (a) The MF-DFA2 fluctuation functions F q (s) are shown 
versus the scale s in a log-log plot for the dyadic random cascade 
model with log-Poisson distribution with parameters A = 1, 5 = 0.9, 
and 7 = — <5 4 /2. The symbols are the same as for Fig. 1. The 
dashed straight lines have the theoretical slopes h(— 10) = 0.743 
and h(+10) — 0.567 and are shown for comparison, (b) The q 
dependence of the generalized Hurst exponent h(q) determined by 
fits is shown for MF-DFA2 and different parameters (see legend). 
The lines are the theoretical values obtained from Eq. (^). In part 
(c) r(q) = qh(q) — 1 is shown. 
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FIG. 4. (a) The modified and rescaled MF-DFA3 fluctuation func- 
tions F q (s)/s[= F q (s)] are shown versus the scale s in a log-log 
plot for a series of independent numbers with a power-law probabil- 
ity density distribution P(x) ~ x~^ a+1 ^ with a = 1. The symbols 
are the same as for Fig. 1. The straight dashed lines have the cor- 
responding theoretical slopes h(— 10) = 1 and h(+!0) = 0.1 and 
are shown for comparison, (b) The q dependence of the generalized 
Hurst exponent h(q) = h(q) — 1 determined by fits on large scales s 
is shown for MF-DFA3 and a = 0.5 (□), 1.0 (O), and 2.0 (A). The 
lines are the theoretical values obtained from Eq. (Q). In part (c) the 
corresponding r(q) is shown. The broad distribution of the values 
leads to multifractality (bi-fractality) in all three cases. 
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FIG. 5. (a) The MF-DFA fluctuation functions Fl h ${s) (open 
symbols) and F^ ut (s) (filled symbols) are shown versus the scale s 
in a log-log plot for randomly shuffled series of long-range correlated 
series with H — 0.75 (□), for the dyadic random cascade model with 
log-Poisson distribution with parameters A = In 1, 5 — 0.9, and 
7 = —6 /2 (O), and for power-law distributed random numbers Xk 
with P(x) ~ x~ 2 (A). The correlations and the multifractality are 
destroyed by the shuffling procedure for the first two series, but for 
the broadly distributed random numbers the multifractality remains. 
The dashed line has the slope H — 0.5 and is shown for comparison, 
(b) The ratios of the MF-DFA2 fluctuation functions F q (s) of the 
original series and the MF-DFA2 fluctuation functions Fq hul (s) of 
the randomly shuffled series are shown versus s for the same mod- 
els as in (a), correlated series (v), dyadic random cascade model 
(O), and power-law distributed random numbers (+ for q = —10, 
x for q — +10). The deviations from the slope h COI — indicate 
long-range correlations. 
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FIG. 6. (a) The q dependence of the slopes /i s h u f(?) of the same 
models as in Fig. |5|(a). The lines indicate the theoretical values: 
H = 0.5 for shuffled data with narrow distribution, and h(q) from 
Eq. ( |27j ) for the series of numbers with a power-law probability den- 
sity distribution. The symbols are the same as in Fig. ^ Part (b) 
shows T a huf (g) = qh Bhui (q) - 1. 
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FIG. 7. (a) Histograms of the generalized Hurst exponents 
h(— 10) (black bars) and ft s huf(— 10) (grey bars) for 100 gener- 
ated monofractal series with H = 0.75. The exponents have 
been fitted to MF-DFA2 fluctuation functions in the scaling range 
400 < s < 2000. (b) Same as (a), but for h(+W) and 
ftshuf (+10). (c,d) Same as (a,b), but for the dyadic random cas- 
cade model with log-Poisson distribution and parameters A = In 1, 
5 — 0.9, and 7 = — <5 4 /2. The corresponding theoretical values are 
h(-W) = 0.743 and h(+W) = 0.567 from Eq.(||) for the original 
series and hshul = 0.5 for the shuffled series. From the histogram 
of h(+10) it would be hard to draw any conclusions regarding mul- 
tifractality. (e,f) Same as (a,b), but for power-law distributed random 
numbers with the distribution P(x) ~ x~ 2 . The corresponding the- 
oretical values from Eq. ( p7j ) are h(— 10) = 1 and h(+10) = 0.1 
for the original and the shuffled series. The length of all series is 
L — 8192. The figure shows that correlations and multifractality due 
to correlations (a-d) are eliminated by the shuffling procedure, while 
multifractality due to a broad distribution (e,f) remains. It further al- 
lows to estimate the statistical fluctuations in the scaling exponents 
h(q) determined by the MF-DFA for monofractal (a,b), correlation 
multifractal (c,d) and distribution multifractal (e,f) series. 
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FIG. 8. (a) The MF-DFA2 fluctuation functions F q (s) are shown 
versus the scale s in log-log plots for the binomial multifractal model 
with a = 0.75 and an additional linear trend Xk — > Xk + k/500L. (b) 
The scaled WTMM partition functions [sZ(q, s)] 1 ^ are shown for 
the same series and the same values of q. The symbols are the same 
as for Fig. 1. (c) The q dependence of the generalized Hurst expo- 
nent h(q) for the generated series with linear trend for the MF-DFA2 
(□) and the second order WTMM (O) methods. Corresponding re- 
sults for a binomial multifractal with an additional quadratic trend 
are also included for MF-DFA2 (A) and second order WTMM (v) 
methods. The quadratic trend causes deviations from the line indicat- 
ing the theoretical values [obtained from Eq. (P0|)], which disappear 
if MF-DFA3 is employed (x). The values of h(q) have been deter- 
mined by fits to F q (s) and Z(q, s) in the regime 50 < s < 2000. 
The relation h(q) = [f(q) + l]/q from Eq. ( |l3| ) has been used to 
convert the exponent f(q) from Eq. (^Tj) into h(q). (d) The q depen- 
dence of r(g). 
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FIG. 10. Same as Fig. g for the binomial model with a = 0.75. 
The theoretical values of the generalized Hurst exponents are 
/i(-10) = 1.90 and ft(+10) = 0.515 according to Eq. @. 
The figure shows that our findings regarding the performance of the 
MF-DFA and WTMM methods for uncorrelated monofractal series 
in Fig. ^ also hold for multifractal series. 



FIG. 9. (a) Histograms of the generalized Hurst exponents 
h(— 10) for 100 random uncorrelated series with H — 0.5. The ex- 
ponents have been fitted to MF-DFA2 fluctuation functions F_io(s) 
(black bars) in the scaling range 40 < s < 2000 and to WTMM 
results Z(—10, s) (grey bars) in the scaling range 5 < s < 250. The 
length of the series is L — 8192. The relation h(q) = [f(g) + l]/q 
from Eq. ( |l3| ) has been used to convert the exponent r(q) from 
Eq. Q into h(q). (b) Same as (a), but for h(+W). (c,d) Same 
as (a,b), but for longer series (L = 65536), where statistical fluc- 
tuations are reduced. The figure shows that the MF-DFA seems to 
give slightly more reliable results than the WTMM method for short 
series and negative moments (q — —10), see (a). In the other cases, 
the performance of both methods is similar. 
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